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Abstract. In this talk we introduce our recent results of global ID MHD simulations for the ac- 
celeration of solar and stellar winds. We impose transverse photospheric motions corresponding 
to the granulations, which generate outgoing Alfven waves. The Alfven waves effectively dissi- 
pate by 3-wave coupling and direct mode conversion to compressive waves in density-stratified 
atmosphere. We show that the coronal heating and the solar wind acceleration in the open mag- 
netic field regions are natural consequence of the footpoint fluctuations of the magnetic fields at 
the surface (photosphere). We also discuss winds from red giant stars driven by Alfven waves, 
focusing on different aspects from the solar wind. We show that red giants wind are highly 
structured with intermittent magnetized hot bubbles embedded in cool chromospheric material. 
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1. Introduction 

The Alfven wave, generated by the granulations or other surface activities, is a promis- 
ing candidate operating in the heating and acceleration of solar winds from coronal holes. 
It can travel a long distance so that the dissipation plays a role in the heating of the solar 
wind plasma as well as the lower coronal plasma, in contrast to other processes, such 
as magnetic reconnection events and compressive waves, the heating of which probably 
concentrates at lower altitude. 

While high-frequency ioncyclotron waves are recently highlighted for the preferential 
heating of minor heavy ions (Axford & McKenzie 1997) the protons, which compose the 
main part of the plasma, are supposed to be mainly heated by low- frequency (< 0.1Hz) 
components in the MHD regime, because (1) the low-frequency wave is expected to have 
more power, and (2) the resonance frequency of the proton is higher than those of heavier 
ions so that the energy of the ioncyclotron wave is in advance absorbed by heavy ions 
(Cranmer 2000). In this paper, we focus on roles of such low frequency Alfven waves. 

When considering low-frequency Alfven waves in solar and stellar atmospheres, the 
stratification of the density due to the gravity is quite important because the variation 
scale of the density, and accordingly Alfven speed, is comparable to or shorter than the 
wavelengths; the WKB approximation is no longer applicable. Also, the amplitude is 
amplified because of the decrease of the density so that waves easily become nonlinear. 
Recently, we have extensively studied the heating and acceleration of solar and stellar 
winds by self-consistent MHD simulations from the photosphere to sufficiently outer 
region (Suzuki & Inutsuka 2005; 2006; hereafter SI06; Suzuki 2007). We review these 
works in this contribution talk. 
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2. Simulation 

We consider ID open flux tubes which are super-radially open, measured by heliocentric 
distance, r. The simulation regions are from the photosphere (density, p — 10 _7 g cm~ 3 , 
for the Sun) to several stellar radii. Radial field strength, B r , is given by conservation 
of magnetic flux as B r r 2 f(r) = const., where f(r) is a super-radial expansion factor(see 
SI06 for detail). 

We input the transverse fluctuations of the field line by the granulations at the photo- 
sphere, which excite Alfven waves. In this paper we only show results of linearly polarized 
perturbations with power spectrum proportional to 1/v, where v is frequency (for circu- 
larly polarized fluctuations with different spectra, see SI06). Amplitude, (dv± t o}, at the 
photosphere is chosen to be compatible with the observed photospheric velocity ampli- 
tude ~ 1km s _1 (Holweger 1978). At the outer boundaries, non-reflecting condition is 
imposed for all the MHD waves, which enables us to carry out simulations for a long 
time until quasi-steady state solutions are obtained without unphysical wave reflection. 

We dynamically treat the propagation and dissipation of the waves and the heating and 
acceleration of the plasma by solving ideal MHD equations. In the energy equation we 
take into account radiative cooling and Spitzer thermal conduction (SI06). We adopt the 
second-order MHD-Godunov-MOCCT scheme (Sano & Inutsuka 2008 in preparation) to 
update the physical quantities. We initially set static atmosphere with a temperature 
T = 10 4 K to see whether the atmosphere is heated up to coronal temperature and 
accelerated to accomplish the transonic flow. At t = we start the inject of the transverse 
fluctuations from the photosphere and continue the simulations until the quasi-steady 
states are achieved. 



3. Results 

3.1. Fast Solar Wind 

Figure 1 plots the initial condition (dashed lines) and the results after the quasi-steady 
state condition is achieved at t = 2573 minutes (solid lines), compared with recent 
observations of fast solar winds. We set the transverse fluctuation, (dv± t o) = 0.7km s _1 , 
and field strength, B r fl = 161G at the photosphere, and the total superradial expansion 
factor, /tot = 75. Figure 1 shows that the initially cool and static atmosphere is effectively 
heated and accelerated by the dissipation of the Alfven waves. The sharp transition region 
which divides the cool chromosphere with T ~ 10 4 K and the hot corona with T ~ 10 6 K 
is formed owing to a thermally unstable region around T ~ 10 5 K in the radiative cooling 
function (Landini & Monsignori-Fossi 1990). The hot corona streams out as the transonic 
solar wind. The simulation naturally explains the observed trend quite well. 

The heating and acceleration of the solar wind plasma in inner heliosphcrc is done by 
the dissipation of Alfven waves. Here we inspect waves in more detail. Figure 2 presents 
contours of amplitude of v r , p, v±, and B±/B r in Rq ^ r 15R Q from t = 2570 min. 
to 2600 min. Red (blue) shaded regions denote positive (negative) amplitude. Above 
the panels, we indicate the directions of the local 5 characteristics, two Alfven , two 
slow, and one entropy waves at the respective positions. Note that the fast MHD and 
Alfven modes degenerate in our case (wave vector and underlying magnetic field are in 
the same direction), so we simply call the innermost and outermost waves Alfven modes. 
In our simple ID geometry, v r and p trace the slow modes which have longitudinal wave 
components, while v± and B± trace the Alfven modes which are transverse 

One can clearly see the Alfven waves in v± and B±/B r diagrams, which have the same 
slopes with the Alfven characteristics shown above. One can also find the incoming modes 
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Figure 1. Results of fast solar wind mode with observations in polar regions. Outflow 
speed, « r (km s _1 ) (top-left), temperature, T(K) (top-right), density in logarithmic scale, 
log(p(g cm~ 3 )) (bottom-right), and rms transverse amplitude, (dv±)(km s" 1 ) (bottom-left) are 
plotted. Observational data in the third panel are electron density, log(iV e (cm~ 3 )) which is to 
be referred to the right axis. Dashed lines indicate the initial conditions and solid lines are the 
results at t = 2573 minutes. In the bottom panel, the initial value ({dv±) = 0) dose not appear. 
The observational data in the inner region (< 67?©) are from SOHO (Teriaca et al.2003; Zangrilli 
et al.2002; Fludra et al.1999; Wilhelm et al.1998; Lamy et al.1997; Banergee et al.1998; Esser et 
al.1999) and those in the outer region are from interplanetary scintillation measurements (Grail 
et al.1996; Habbal et al.1994; Kojima et al.2004; Canals et al.2002). 



propagating from lower-right to upper-left as well as the outgoing modes generated from 
the surface These incoming waves arc generated by the reflection at the 'density mirrors' 
of the slow modes. At intersection points of the outgoing and incoming characteristics the 
non-linear wave-wave interactions take place, which play a role in the wave dissipation. 

The slow modes are seen in v r and p diagrams. Although it might be difficult to 
distinguish, the most of the patterns are due to the outgoing slow modest which are 
generated from the perturbations of the Alfven wave pressure, B±/8ir (Kudoh & Shibata 
1998& Tsurutani et al. 2002). These slow waves steepen eventually and lead to the shock 
dissipation. 

Figure 3 presents the dissipation of the waves more quantitatively. It plots the following 
quantities, 

of outgoing Alfven , incoming Alfven , and outgoing slow MHD (sound) waves, where 

f The phase correlation of the longitudinal slow waves is opposite to that of the transverse 
Alfven waves. The outgoing slow modes have the positive correlation between amplitudes of v r 
and p, (8v r Sp > 0), while the incoming modes have the negative correlation (5v r 8p < 0). 
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Figure 2. r — t diagrams for v r (upper-left), p (lower-left), v± (upper-right), and B±/B r (low- 
er-right.) The horizontal axises cover from Rq to 15-R©, and the vertical axises cover from 
t = 2570 minutes to 2600 minutes. Dark(blue) and light shaded regions indicate positive and 
negative amplitudes which exceed certain thresholds. The thresholds are dv r — ±96km/s for v r , 
dp/ p — ±0.25 for p, v± = ±180km/s for v±, and B±/B r — ±0.16 for B±/B r , where dp and dv r 
are differences from the averaged p and v r - Arrows on the top panels indicate characteristics of 
Alfven , slow MHD and entropy waves at the respective locations. 

Sv and v p h are amplitude and phase speed of each wave mode. S c is an adiabatic con- 
stant derived from wave action (Jacques 1977) in unit of energy flux. For the incoming 
Alfven wave, we plot the opposite sign of S c so that it becomes positive in the sub- 
Alfvenic region. The outgoing and incoming Alfven waves are decomposed by correlation 
between v± and B± . Extraction of the slow wave is also from fluctuating components of 
v r and p. 

Figure 3 clearly illustrates that the outgoing Alfven waves dissipate quite effectively; 
S c becomes only ~ 10~ 3 of the initial value at the outer boundary. A sizable amount 
is reflected back downward below the coronal base (r — i?s < 0.01i?s), which is known 
from the incoming Alfven wave following the outgoing component with slightly smaller 
level. This is because the wave shape is considerably deformed owing to the steep density 
gradient; a typical variation scale (< 10 5 km) of the Alfven speed becomes comparable 
to or even shorter than the wavelength (= 10 4 — 10 6 km). Although the energy flux, 
~ 5 x 10 5 erg cm~ 2 s -1 , of the outgoing Alfven waves [S c in the static region is equivalent 
with the energy flux) which penetrates into the corona is only ~ 15% of the input value, 
it satisfies the requirement for the energy budget in the coronal holes (Withbroe & Noyes 
1977). 

The processes discussed here are the combination of the direct mode conversion to 
the compressive waves and the parametric decay instability due to three-wave (outgoing 
Alfven , incoming Alfven , and outgoing slow waves) interactions (Goldstein 1978; Tera- 
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Figure 3. S c of outgoing Alfven mode (solid), incoming Alfven mode (dotted), and outgoing 
MHD slow mode (dashed) at t = 2573 mins. Hatched region indicates the chromosphere and 
low transition region with T < 4 x 10 4 K. 

sawa ct al. 1986) of the Alfven waves. These processes, which are not generally efficient in 
homogeneous background, become effective by amplification of velocity amplitude in the 
density decreasing atmosphere. The Alfven speed also varies a lot even within one wave- 
length of Alfven waves with periods of minutes. This leads to both variation of the wave 
pressure in one wavelength and partial reflection through the deformation of the wave 
shape (Moore et al.1991). The density stratification plays a key role in the propagation 
and dissipation of the Alfven waves. 

3.2. Evolution to Red Giant Winds 

So far we have focused on the acceleration of solar wind. The same process is expected to 
operate in other types of stars that have surface convective layer, such as red giant stars, 
proto-stars, and intermediate and low mass main sequence stars. As a demonstration, we 
have applied our solar wind simulations to red giant winds (Suzuki 2007). 

We consider stellar winds from 1M@ stars in various evolutionary stages from main 
sequence to red giant branch. The properties of surface fluctuations (e.g. amplitude and 
spectrum) can be estimated from conditions of surface convection which depend on sur- 
face gravity and temperature (e.g. Renzini et al.1977; Stein et al. 2004). Then, we carry 
out the simulations of the red giant winds in a similar manner to the solar wind simula- 
tions. 

Figure 4 presents the evolution of stellar winds of a 1M star from main sequence to 
red giant stages. The middle panel shows that the average temperature drops suddenly 
from T ~ 7 x 10 5 K in the sub-giant star (blue) to T ^ 10 5 K in the red giant stars, which 
is consistent with the observed "dividing linc"(Linsky & Haisch 1979). The main reason 
of the disappearance of the steady hot coronae is that the sound speed (« 150 km s _1 ) 
of w 10 6 K plasma exceeds the escape speed, v csc {r) = y / 2GM±/r, at r > a few i?* in 
the red giant stars; the hot corona cannot be confined by the gravity any more in the 
atmospheres of the red giant stars. Therefore, the material flows out before heated up to 
coronal temperature. 

In addition, the thermal instability of the radiative cooling function (Landini & Monsignori- 
Fossi 1990) plays a role in the sudden decrease of temperature. Because of the thermal 
instability, magnetized hot (> 10 6 K) bubbles intermittently exist in red giant winds, 
while most of the wind material consist of cool (< 10 4 K) chromospheric gas (Suzuki 
2007). 
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Figure 4. Time- averaged stellar wind structure of the IMq stars. From the top to the bottom, 
radial outflow velocity, v r (km s _1 ), temperature, T (K), and density, p(g cm -3 ), are plotted. 
The black, blue, green, and red lines are the results of stellar radii, R = Rq (the present Sun), 
3.17?© (sub-giant), WRq (red giant), and 317?© (red giant), respectively. 



4. conclusions 

We have performed ID MHD numerical simulations of solar and stellar winds from the 
photosphere. The low-frequency Alfven waves are generated by the footpoint fluctuations 
of the magnetic field lines. We have treated the wave propagation and dissipation, and 
the heating and acceleration of the plasma in a self-consistent manner. Our simulation 
is the first simulation which treats the wind from the real surface (photosphere) to the 
(inner) heliosphere with the relevant physical processes. 

We have shown that the dissipation of the low-frequency Alfven waves through the 
generation of the compressive waves (decay instability) and shocks (nonlinear steepen- 
ing) is one of the solutions for the heating and acceleration of the plasma in the coronal 
holes. However, we should cautiously examine the validity of the 1-D MHD approximation 
we have adopted. There are other dissipation mechanisms due to the multidimensional- 
ity, such as turbulent cascade into the transverse direction (Goldrcich & Sridhar 1995; 
Oughton ct al.2001) and phase mixing (Heyvaerts & Priest 1983). If Alfven waves cascade 
to higher frequency, kinetic effects (e.g. Nariyuki & Hada 2006) becomes important. 

We have also extended the solar wind simulations to red giant winds. With stellar evo- 
lution, the steady hot corona with temperature, T w 10 6 K, suddenly disappears because 
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the surface gravity becomes small; hot plasma cannot be confined by the gravity. Ther- 
mal instability also generate intermittent magnetized hot bubbles in cool chromospheric 
winds. 
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